Magneto-optical Trapping of Cadmium 
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We report the laser-cooling and confinement of Cd atoms in a magneto-optical trap, and charac- 
terize the loading process from the background Cd vapor. The trapping laser drives the 1 5o — » 1 P\ 
transition at 229 nm in this two-electron atom and also photoionizes atoms directly from the 1 Pi 
state. This photoionization overwhelms the other loss mechanisms and allows a direct measurement 
of the photoionization cross section, which we measure to be 2(1) x 10~ 16 cm 2 from the 1 Pi state. 
When combined with nearby laser-cooled and trapped Cd + ions, this apparatus could facilitate 
studies in ultracold interactions between atoms and ions. 
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INTRODUCTION 



The magneto-optical trap (MOT) is an indispensable 
source of cold atoms for a range of studies and appli- 
cations in atomic physics, from precision atomic spec- 
troscopy [H and cold collisions Q to atom interferome- 
try and the generation of quantum-degenerate gases 
While nearly all cold atom experiments deal with the al- 
kali atoms, there has been progress in the trapping of 
two-electron atomic species such as Ca, Mg, Sr, and Yb 
0, 0, S 0)j mainly for experiments involving high res- 
olution spectroscopy of the 1 S — > 3 P intercombination 
lines. We report here the trapping of atomic Cd atoms 
in a deep- ultraviolet MOT operating on the 1 Sq — > 1 P\ 
transition at 229 nm. 

When producing a Cd MOT, the trapping light can 
also photoionize the atoms directly from the 1 P\ excited 
state. While this introduces losses on the trapping pro- 
cess, it also provides an opportunity to reliably create 
cold ions and atoms at the same location 0] for the in- 
vestigation of ultracold atom- ion interactions @,[l0]. One 
interesting future possibility is the transfer of coherence 
between ground state hyperfine levels in a trapped ion to 
pure nuclear spin states in a neutral atom lacking elec- 
tron spin. Because the nuclear spin can be extremely 
well-isolated from environmental influences 111 ], control 
of such a coherent transfer process may have applications 
to the long term storage of quantum information. 

In this work we realize the first Cd MOT and charac- 
terize the various trapping parameters. Results are com- 
pared with simple analytic and Monte-Carlo simulation 
models of the trapping process. Through a detailed in- 
vestigation of the loss rate as a function of laser intensity, 
the absolute photoionization cross section from the 1 P\ 
state is determined. 



BACKGROUND 

Cadmium has eight stable isotopes, six of which are 
relatively abundant. Fig. [1] shows the electronic struc- 
ture of Cd for both bosons (nuclear spin 1=0, even iso- 
topes) and fermions (1=1/2, odd isotopes). Most of the 
the data presented here is for 112 Cd. The 1 Sq - 1 P\ 
atomic transition used for the MOT occurs at a wave- 
length of A=228.8 nm with an excited state lifetime of 
r=1.8 ns (radiative linewidth 7/27r = 91 MHz) and sat- 
uration intensity of I sa t—^hc^/(3\ 3 ) ~ 1.0 W/cm 2 . The 
saturated photon recoil acceleration on a Cd atom is 
a = h/2Tm\ = 4.4 x I0 5 g, which is 50 times that of Rb 
(here g is the acceleration due to gravity and m is the 
mass of a single Cd atom). Note that the 228.8 nm light 
can also excite atoms from the 1 Pi state directly to the 
ionization continuum. 

In a vapor cell, the radiative forces accumulate atoms 
following the rate equation 



(i) 



where N is the number of trapped atoms, L is the loading 
rate, T is the loss rate related to single atom effects, (3 is 
the loss rate due to binary collisions within the trap, and 
V is the effective volume occupied by the trapped atoms 

B 



13, 14, HI. Using simple kinetic gas theory at con- 
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stant temperature one can show that L « nV c /v 
where V c is the capture volume, v c is the capture velocity 
[13], vth is the thermal velocity, and n is the density of 
Cd atoms in the background vapor 0, [ijj ■ For com- 
parison to the data we use a simple analytic 1-D laser 
cooling model to find the capture velocity, as detailed 
in appendix 1. Appendix 2 presents a 3-D Monte-Carlo 
simulation, which includes magnetic field and polariza- 
tion effects, to directly estimate the loading rate. 

When the MOT density is low (< 10 9 atoms/cm 3 ), 
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FIG. 1: Cadmium energy level diagram (not to scale), a. 
The bosonic (even) isotopes (1=0) of Cd. b. The fermionic 
(odd) isotopes (1=1/2) of Cd, where the 1 P\ hyperfine split- 
ting arises from (L ■ I) coupling. Individual levels are labeled 
with mp. 



the atoms are essentially non-interacting and we expect 
the density to be limited by the cloud temperature. In 
this regime the spatial distribution of trapped atoms is 
expected to be Gaussian with a cloud radius that is in- 
dependent of the trapped atom number. This contrasts 
with high density (>10 10 atoms/cm 3 ) MOTs where ef- 
fects such as reradiation [l8| must be considered. The Cd 
MOT reported here operates in the low density regime, 
and the last term of Eq. [TJcan be neglected. Unlike con- 
ventional alkali MOTs, where single atom loss mecha- 
nisms primarily involve collisions between trapped atoms 
and the background gas, Cd (like Mg) has an additional 
single atom loss term due to photoionization [3, l20j ]. 
Solving Eq. [1] for the steady state number of trapped 
atoms gives N ss = L/T, with the loss rate given by 



r = r + r. 



(2) 



Here Tq represents the rate at which trapped atoms 
are ejected due to collisions with the background vapor 
(dominated by Cd) and r.; on is the photoionization rate: 



r. 



aP(I, 8)1 



(3) 



In this expression, a is the photoionization cross section, 
fku is the photon energy, / is the total MOT laser beam 
intensity, and P(I, S) is the fraction of atoms in the ex- 
cited state ( 1 Pi) defined as 



P(I,S) 
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where S = A/7 is the laser detuning scaled to the natural 
linewidth and s = I j I sa t is the saturation parameter. 
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FIG. 2: Left: Schematic diagram of the laser system and 
the laser lock (DAVLL). The laser lock consists of the Cd 
cell, a quarter wave plate (QWP), a polarizing beam splitter 
(PBS), and two photodetectors (DA, DB) for path a and b, 
respectively. Right: The MOT vacuum chamber and the laser 
beam geometry. The MOT is formed by 6 independent beams. 
The imaging system sits below the chamber, and the dark 
shaded regions are the NdFeB magnets. 



EXPERIMENTAL SET-UP AND PROCEDURE 

A schematic of the experimental apparatus is shown 
in Fig. O Since Cd has a large linewidth, high magnetic 
field gradients are required to shift the Zccman levels suf- 
ficiently for the atoms to feel a substantial trapping force 
at the edge of the laser beams. We use NdFeB permanent 
ring magnets with a 2.54 cm outer diameter, 0.64 cm in- 
ner diameter, and 0.95 cm thickness that are mounted 
coaxially on translational stages. By adjusting the axial 
separation of the magnets we can achieve magnetic field 
gradients up to 1500 G/cm at the trap center. 

The trapping beams are generated with a frequency 
quadrupled ThSapphire laser, yielding 2.5 mW at 228.8 
nm. The ultraviolet light is split into six independent 
trapping beams in order to better control the intensity 
balance of the counter-propagating beams. The MOT 
can withstand an intensity imbalance of 10% between a 
pair of beams, and we can balance the intensity between 
any pair of counter-propagating beams to better than 
5%. Typical beam waists range from w =0.5 mm to 
1.5 mm and the total power ranges from P= 0.7 mW to 
2.0 mW, resulting in peak intensities ranging from about 
0.03 W/cm 2 to 0.5 W/cm 2 . 

Approximately 200 /iW is split from the main laser 
beam and directed to a small cadmium vapor cell to sta- 
bilize the laser frequency. We use a dichroic atomic va- 
por laser lock (DAVLL) [HI, ^ operating on the l S 
— > 1 P\ transition in Cd. The cell is heated to 80° C to 



3 



increase optical absorption to about 80% through the 5 
cm cell. A uniform magnetic field is applied along the 
laser beam axis to lift the degeneracy of the 1 P 1 states. 
When linearly polarized light is sent through the cell the 
difference between absorption of the Zeeman-shifted a + 
and er~ transitions produces a dispersive-shaped signal 
and the laser is locked to the zero crossing point of this 
signal. The capture range is determined by the Zceman 
splitting between the two transitions, or about 1.5 GHz 
in a 500 G field. To change the detuning we move the zero 
crossing point by attenuating the laser power in one of 
the polarization paths (a or b) after the cell (see Fig. [2]). 
The lock is stable to within 30 MHz, or 0.37, over the 1.5 
GHz capture range, and the dominant sources of fluctu- 
ations are beam-steering drifts and birefringence fluctu- 
ations of the cell windows from temperature drifts over 
times greater than 1 second. 

The MOT chamber contains a 1 cm long hollow stain- 
less steel tube of diameter 0.1 cm packed with about 0.02 
g of pure Cd wire. We control the background Cd vapor 
pressure throughout the entire chamber by heating this 
small oven. When we direct the trapping beams into 
the chamber we see tracks of fluorescing Cd within the 
extent of the laser beams. Based on this atomic fluores- 
cence, we estimate the background Cd vapor pressure to 
range between approximately 10~ n torr with the oven 
off to about 10~ 10 torr with the oven at approximately 
300° C. We speculate that the Cd atoms sublimated from 
the oven do not readily stick to the chamber surface, re- 
sulting in good control of the Cd vapor pressure with the 
small oven. We note that the vapor pressure of Cd is pre- 
dicted to be 10~ n torr at room temperature [Hj], which 
is consistent with our observations. 

The atomic fluorescence from the trapped atoms is col- 
lected with an f/3 lens (a solid angle of dil/4n — 0.6%) 
and imaged onto an intensified charge coupled device 
(ICCD) camera. Every photon incident on the cam- 
era yields 77G ~ 65 counts, where 77=20% is the quan- 
tum efficiency of the camera and G is the ICCD gain 
factor. Including an optical transmission of T ~ 50% 
in the imaging system, we expect a total count rate 
of jP(I,5)GriT(d£l/4:iT) ~ 10 7 counts/sec from each 
trapped atom in the MOT. In this way, we relate the to- 
tal fluorescence count rate to the number of atoms in the 
MOT, with an estimated accuracy of 50%. For various 
settings of the MOT parameters, we are able to observe 
between about 10-3000 atoms in the MOT. 



RESULTS AND DISCUSSION 

A typical observation of the fluorescence growth from 
trapped atoms in time is shown in Fig. 02 allowing a 
determination of the steady-state number of atoms and 
the net loss rate, T, from the trap. An image of the 
fluorescence distribution from the trapped atoms is also 
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FIG. 3: Left: Typical loading curve showing the buildup 
in the MOT fluorescence vs. time. For this data set, the 
MOT parameters are laser power P=1.45 mW, beam waist 
w=1.25 mm, detuning S=— 0.7, and magnetic field gradient 
B'=500 G/cm. The steady state MOT number is calculated 
from the fluorescence signal and for this data the buildup 
time is 1.5 sec. Right: MOT image taken with the camera for 
N ss = 1200 atoms. The MOT parameters for this data set are 
P =1.45 mW, w=2.5 mm, <5=— 0.7, and B'=500 G/cm. The 
integration time for the camera was 5 ms. A 2-D Gaussian fit 
to the image yields an rms radius of 200 )im and a peak atom 
density of 10 s atoms/cm 3 . 



shown in Fig.[3l revealing a Gaussian-shaped atom cloud 
as expected from the temperature-limited density. The 
typical geometric mean rms radius of the MOT is 200 /im, 
with some dependence upon the magnetic field gradient, 
laser power and detuning. The largest MOT we have 
observed held approximately 3000 atoms, with a peak 
density of about 10 s atoms/cm 3 . 

Fig. 2] shows the steady state number of atoms, N ss , 
in the MOT vs. magnetic field gradient, B' , for beam 
waist u>=1.25 mm, detuning 5=— 0.6, and a total power 
P=1.8 mW. Under these conditions the maximum steady 
state number is observed at 500 G/cm. At this optimal 
field gradient, the Zeeman shift of the excited state lev- 
els at the edge of the laser beam is approximately one 
linewidth. Above this optimal value the steep magnetic 
field gradient shifts the atoms out of resonance with the 
laser beams, reducing the capture volume. At lower field 
gradients N ss quickly decreases, presumably due to a 
lower trap depth resulting from an increased sensitivity 
to trapping parameters. 

From the equipartition theorem we obtain a relation 
connecting the cloud radius and temperature, nr 2 = k(,T, 
where r is the atomic cloud rms radius, kb is Boltzmann's 
constant, T is the temperature in Kelvin, and k is the 
trap spring constant n — 8fibksB'5/(l + s + 4<5 2 ) [2^ ]. In 
this expression, fib is the Bohr magneton, and k = 2tt/\ 
is the wavenumber. Replacing T with the Doppler tem- 
perature, To = firf(l + s + 4S 2 ) / (8kb\S\) , gives a relation 
between the temperature-limited cloud radius and the 
magnetic field gradient: 
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FIG. 4: Observed steady-state MOT number vs. axial mag- 
netic field gradient B' (points), along with the 3-D model 
(solid line) for P=0.8 mW, 5=— 0.6, and w=2.5 mm. 
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FIG. 6: Observed steady-state atom number vs. S (points) 
along with the 1-D (dotted line) and 3-D (solid line) models 
for P=l.8 mW, B'=500 G/cm and w=2.5 mm. 
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FIG. 5: Atom cloud rms diameter vs. B' for P=0.8 mW, 
S——0.6, and w=2.5 mm. A characterization is provided by 
the longest (circles) and shortest (squares) rms size of the el- 
liptical MOT. The diameter is about 5 times larger than what 
Doppler theory predicts. The solid lines show the (B') -1 / 2 
dependence expected from Eq. [5] 



T ~ y 64^ b S 2 ksB' ' ( ' 

Fig. [5] shows the MOT rms radius vs. magnetic field 
gradient; as expected from Eq. [5l the cloud gets smaller 
as B' increases. The MOT diameter is roughly 5 times 
larger than what Doppler theory predicts. Similar results 
were found in Sr, where the MOT temperature exceeded 
the expected Doppler temperature 0]. 

The dependence of the steady-state number of trapped 
atoms on MOT detuning and laser power is shown in 
Figs. [6] and [7] In both figures, the experimental data 
is plotted along with the 1-D and 3-D theoretical pre- 
dictions. The observed number of trapped atoms is 1- 
2 orders of magnitude below predictions, likely due to 
alignment imperfections and intensity imbalances not in- 
cluded in the models. Fig. [8] shows how the measured 
atom cloud size decreases as the MOT laser power is in- 
creased (at a fixed beam waist), as expected from Eq. [5] 

In Fig. the filling of the MOT is shown for Cd va- 
por pressures of approximately 10 -10 torr and 10 -11 torr. 
Unlike conventional vapor cell MOTs, we find that the 
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FIG. 7: Observed steady-state atom number vs. power 
(points) for 5=— 0.7, £>'=500 G/cm and w=2.5 mm along 
with the 1-D (solid line) and 3-D (dotted line) models. 
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FIG. 8: MOT cloud diameter vs. total MOT laser power for 
S=— 0.6, B'=500 G/cm and w=2.5 mm. The solid lines show 
the expected dependence of the MOT diameter on power from 
Eq.0 



filling time (loss rate) is independent of the background 
pressure, while the steady-state number of atoms in the 
MOT is strongly dependent on pressure. This indicates 
that collisions with the background gas have very little 
effect on the loss rate and instead we are limited by pho- 
toionization loss from the MOT beams. This is investi- 
gated in more detail by measuring the filling time (loss 
rate) as the MOT laser intensity is varied, as shown in 
Fig. [TUJ We observe a roughly quadratic dependence of 
loss rate on intensity, consistent with Eq. [5] The extrap- 
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olated loss rate at zero intensity is much smaller than all 
of the observations, directly indicating that Tq <C r\ on , 
or that the loss rate in this experiment is dominated by 
photoionization. From this measurement, we can also 
directly extract the photoionization cross section from 
the 1 P\ state, given measurements of the intensity, ex- 
cited state fraction P(I,5) , and the known wavelength of 
the light. We find that the photoionization cross sec- 
tion of the 1 Pi state of Cd from the 228.8 nm light is 
a = 2(1) x 1CP 16 cm 2 , with the error dominated by un- 
certainties in the laser intensity and detuning. 
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FIG. 10: Observed loading rate vs the saturation parameter 
s=I/I sa t- The power is varied for a constant beam waist of 
w=1.25 mm. The photoionization cross section out of the 
x Pi state is determined from a quadratic fit to s given by 
Eq. [2] Extrapolating the curve to zero intensity (not shown 
here) gives information on the loss rate due to collisions with 
background gas. 
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FIG. 9: Top: Observed trapped atom number N(t) for two 
different Cd background vapor pressures. The top curve cor- 
responds to a pressure of 10~ 10 torr and the lower curve cor- 
responds to 1CP 11 torr. By fitting the data to a growing ex- 
ponential, N(t) = iV ss (l — e _rt ), we find that the filling time, 
F _1 , is approximately 1 sec for each case. This is clear from 
the lower logarithmic plot of the data. Bottom: N ss -N(t) 
plotted for both vapor pressures on a log scale. The filling 
times are about 1 sec for each curve. 



FIG. 11: Top: Scan across frequency showing the different 
Cd isotope MOTs. The underlying curve is the Doppler fluo- 
rescence profile of the Cd atoms. At certain frequencies there 
is a large build up, due to the MOT accumulation as its res- 
onance is crossed. Bottom: Natural abundance of neutral 
cadmium isotopes. Out of these eight isotopes, we are only 
able to clearly observe trapping of the four most abundant 
bosonic (even) isotopes. 



FERMIONIC ISOTOPES 

Scanning the laser frequency allows cooling and trap- 
ping of different cadmium isotopes, as shown in Fig. 1111 
We observe that the peak heights correspond to the nat- 
ural abundance of each isotope, showing that the bosonic 
isotopes are equally capable of being trapped. However, 
there is a lack of evidence for the fcrmionic isotopes being 
loaded in the MOT. This is due to the hyperfine struc- 
ture present in the 1 P± states of the fermionic isotopes. 
As shown in Fig. QJ>, the two excited hyperfine states 
for both m Cd and 113 Cd are separated by about 300 
MHz, which is comparable to the natural linewidth of 
Cd. A laser tuned to the red of the upper hyperfine state 



(F'=3/2) but to the blue of the lower hyperfine state 
(F'=l/2) may drive excessive transitions to the lower ex- 
cited state, which could result in too much heating and 
prevent trapping. In addition, the optical transitions be- 
tween the F=l/2 ground states and F'=l/2 excited states 
do not result in spatially dependent differential optical 
pumping by a + and cr~ transitions, a necessary condi- 
tion for a standard MOT. Similar results were reported 
for Yb Q, where much smaller or no MOT was observed 
for fermionic isotopes. In the present case it could be 
that there is a very small fermionic MOT being formed 
but it is not resolvable from the background noise. It may 
be possible to laser cool and trap fermionic isotopes with 
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a dichroic MOT [25(. Here, the cooling laser is tuned to 
the red of the lower hyperfine transitions (F'=l/2) to pro- 
vide the major scattering force for laser cooling and then 
a small fraction of the laser power is frequency shifted to 
the red of the upper hyperfine state (blue of the lower 
state). When this second laser beam is collimated with 
a smaller beam waist, and overlapped with the beam of 
the first color, the laser cooled atoms can be trapped 
in the MOT by driving the upper transitions (F=l/2 to 
F'=3/2). Alternatively, one can work in a much higher 
magnetic field gradient to overwhelm the excited state 
hyperfine structure. In this Paschen-Bach regime, one 
will drive J = to J = 1 transitions to produce a MOT. 
Given a beam waist of 1.0 mm, the required field gradi- 
ent for the MOT will be on the order of 10 4 G/cm, which 
can be realized by a pair of needle electromagnets [26j |. 
The capture volume of the MOT will be much smaller, 
but this scheme may still be useful for single-atom MOT 
experiments. Another alternative is to use a higher laser 
power allowing one to tune to the red of both hyperfine 
states. With a larger detuning, \8\ ^S> Shf, the optical ex- 
citation to the lower and upper manifolds is driven more 
evenly and can produce both cooling and trapping forces 
for the atoms. 



CONCLUSION 



APPENDIX 

1-D DERIVATION FOR STEADY-STATE 
NUMBER OF ATOMS COOLED TO REST IN A 
VAPOR CELL 

The following appendix estimates the number of atoms 
cooled to rest in a vapor cell. For simplicity we assume 
the laser beams to have a top-hat profile. 

The force on an atom moving with velocity, v, in two 
counter-propagating laser beams is 



1 + s + 4(5 - u) 2 



1 



•4(5- 



(6) 

where the scaled velocity is defined as u^kv/j. 

To find the capture velocity v c , that is the maximum 
velocity an atom can possess and still be slowed to rest 
within the cooling laser beams, we must integrate the 
velocity-dependent acceleration a(v) 



vdv 



XQ+l 



dx. 



(7) 



Here / is the laser beam diameter, x Q is defined to be the 
edge of the laser beam, and 



In this paper we present the first Cd magneto-optical 
trap, operating on the 1 Sq - 1 P\ transition at 228.8 nm. 
A characterization of the MOT as a function of mag- 
netic field gradient, detuning, and intensity is presented. 
The same beams that form the MOT also photoionizc 
the atoms inside the MOT. We observe photoionization 
as the dominant loss mechanism, and characterize the 
photoionization cross section. 

This system, when combined with cold ions, opens the 
possibility of studying ultra-cold charge exchange colli- 
sions. One outcome of these studies is the possible trans- 
fer of coherent information from the ion to the neutral 
atom. A possible experiment is to prepare the ion in a 
quantum superposition of the hyperfine qubit states and 
then allow the ion to undergo an ultracold charge ex- 
change with a nearby neutral atom. This results in the 
charge neutralization of the ion, but could also leave some 
of the previously prepared quantum information intact 
in the nucleus. This could allow quantum information to 
be carried by pure nuclear spins with very little interac- 
tion with the environment. Subsequent coherent charge 
exchange with another ion could then allow the nuclear 
quantum information to be manipulated and processed 
using conventional ion trap techniques. In addition to 
applications for quantum information, the long-lived 3 Po 
state could be of interest for optical clocks [l[ and the nar- 
row linewidth of the 1 S'o- 3 Pi transition (70 kHz) would 
allow for an extremely low cooling limit 



a(v) = — — 



1 + s + 4(5 - u) 2 1 + s + 4(5 + u) 2 
Solving Eq. [7] gives 



(8) 



-16hk 3 sSl 16 5 8., ..», 3 ,.0,2 

2 = ^< + -{\ + s -A5 2 )ul + {l + s + A5 2 ) 2 u c . 

(9) 

This fifth order polynomial can be solved numerically to 
find the capture velocity, v c — ju c /k. 

From the steady state solution to Eq. [T] given above, 
we get 



fnhiol 2 v\ 



v\ h hioG c n + v 3 h a ion IP(I, 5) 



(10) 



where a c is the collision cross section, Ui on is the pho- 
toionization cross section, and / is the relative abundance 
of the isotope of interest. 



3-D DERIVATION FOR STEADY-STATE 
NUMBER IN A VAPOR CELL MOT 

A 3-dimensional computer simulation that includes the 
effects of the magnetic field to obtain the steady-state 
number of atoms trapped in an even-isotope Cd MOT 
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is employed to compare with our data. Previous work 
has used models that calculate the capture velocity for 
an individual atom and extract the loading rate from v c 



271 1281. 1291 . 1301 ] or has included individual photon recoil 
events |31| . 

Here, we treat the atoms as non-interacting, point par- 
ticles and examine the dynamics of an ensemble of indi- 
vidual atoms under the application of the laser radia- 
tion. The atoms are subjected to a time-averaged force, 
i.e. we do not track individual photon absorption and 
re-emission events and instead we calculate the averaged 
momentum kicks over hundreds of scattering events. The 
atoms' motional behavior is found by numerically inte- 
grating the position- and velocity-dependent radiation 
force. 

Specifically, taking the ^-direction as an example, the 
net acceleration of an individual atom is given by the sum 
of the acceleration due to the +£ and — x laser beams: 



a x (r,v x ) = - — s x (r) 
2m 



Y l + s tot (f) + (2S+, q (r,v x )) i 



^ l + s tot (r) + (25„, q (r,v x )) 2 



(11) 



where q is the polarization index, and the ± subscripts 
correspond to the +x- and — i-direction beams, respec- 
tively. 

The individual beam saturation parameter, s x , and the 
total saturation parameter, s t ot, are given by: 



s x (r) = — e *"o 

J sat 



stot(r) = 2s x {f) + 2s y (f) + 2s z {f), 



(12) 
(13) 



where I x is the intensity of the ±x laser beams (here 
assumed to be balanced) and s y (f) and s z (r) are denned 
analogously to s x (r) . 

The fraction of the incoming laser radiation that the 
atom experiences as a^- or 7r-polarized is: 



flh zp 1 xB ' 



q=-l (a-) 

1 - (p±,-i +p±,+i), q = (?r) 



P±, g (r) 



where B(r) = B'y z 2 + \{x 2 + y 2 ) is the magnitude of 
the magnetic field written in terms of the magnetic field 
gradient, B 1 , along the strong (z) axis. 

The effective detuning for the atomic transition is given 
by: 

S±, q (r,v x ) = (AT kv x )h + q ^ 9 ^ (15) 

where gp = 1 is the Lande g- factor. Note that for 
the even isotopes of Cd, there is no Zeeman shift for 
7r-polarized radiation. 



Atoms are initially placed uniformly distributed in po- 
sition within a simulation volume with a lateral dimen- 
sion of 6w (3 beam diameters) and the total number of 
atoms is chosen to correspond to the background den- 
sity of the atomic vapor. The atoms are given initial 
velocities distributed according to a Maxwell-Boltzmann 
velocity distribution centered about 0. To save compu- 
tation time, the atoms are discarded if they have an ini- 
tial velocity Vi > 5v c , where v c w 20m/s is the maxi- 
mum capture velocity for the atoms calculated from the 
1-dimensional model (appendix A). For select parame- 
ter sets, we have checked this time-saving assumption 
against a version that simulates all atoms regardless of 
velocity and found no difference in the obtained results. 

For each time step, dt, the time-averaged acceleration 
on each atom is computed (Eq. [TTj) and the new velocity 
and position are calculated according to: 



a T dt 



x = x + v x dt + — a x (dty 



(16) 
(17) 



At each time step, any atoms that have left the simula- 
tion volume are eliminated. The number of new atoms 
added at each time step is calculated by considering the 
average number of atoms that would leave the volume at 
each time step assuming no radiation forces were present. 
These new atoms are added with uniformly distributed 
positions along one of the (randomly) chosen simulation 
box edges and are given velocities that point inwards but 
correspond to a Maxwell-Boltzmann distribution (again 
discarding velocities greater than 5v c ). 

Every 100 time steps, we count the number of atoms 
with positions that are within wq/2 of the origin. As a 
function of time, this sum produces a linear curve with a 
slope that corresponds to the loading rate of atoms being 
trapped by the MOT. The steady-state number of atoms 
confined by the MOT can then be calculated by includ- 
ing isotope abundance and the background collision and 
photoionization loss rates (similarly to Eq. [IT))) . 
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